Local scale invariance as dynamical space-time symmetry in phase-ordering kinetics 
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The scaling of the spatio-temporal response of coarsening systems is studied through simulations 
of the 2D and 3D Ising model with Glauber dynamics. The scaling functions agree with the pre- 
diction of local scale invariance, extending dynamical scaling to a space-time dynamical symmetry. 
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Ageing phenomena are observed in a broad variety of 
systems with slow relaxation dynamics. Useful insight 
may often be gained via the consideration of simple ferro- 
magnetic models (rather than genuinely glassy systems) , 
quenched to a temperature T < T c from a disordered 
initial state into the low-temperature phase. For their 
study two-time quantities, rather than the more usual 
one-time quantities, are particularly useful, see [1-3] for 
recent reviews. Those display dynamical self-similarity 
in the ageing regime where the order-parameter auto- 
correlator C(t, s) = (cj)(t)(f>(s)) decays from its plateau 
value 5ea = to zero (M eq is the spontaneous mag- 
netization). In this regime both s and r = t — s are 
much larger than the microscopic time scale (set to 
unity). Similarly, the autoresponse function, defined as 
R(t,s) — 5((f>(t))/8h(s) with t > s shows a scaling be- 
haviour such that the scaling laws 



C(t,s)~M*f c (t/s) , R(t,s)~s- 1 - a f R (t/s) (1) 



with /c(l) = 1 are found to hold for a broad range of 
models [1-3]. For x 3> 1, the scaling functions usu- 
ally fall off as f c (x) ~ x~ Xc l z and f R (x) ~ x~ x ^ z , 
where z is the dynamic critical exponent, Ac is the au- 
tocorrelation exponent [4,5] and Xr is the autoresponse 
exponent [6]. We shall focus here on phase-ordering 
in the Glauber-Ising model in d > 1 dimensions where 
a = 1/z = 1/2 is expected [7,2]. Recent arguments [8] 
leading to a = 1/4 in the 2D Glauber-Ising model have 
been rejected through a detailed study of the scaling of 
R(t, s), the results of which reconfirm a — 1/2 [9]. 

Eq. (1) states that the two-time quantities evaluated 
at the same spatial location transform covariantly under 
a global rescaling of time t — > bt with b constant. Re- 
cently, it has been proposed that the response functions 
should transform covariantly under local scale transfor- 
mations with b — b(t), but with time-translations ex- 
cluded [10]. By analogy with conformal invariance, in 
particular covariance of R under the so-called 'special' 
transformations which transform time as t —> t/(l + "ft), 
is assumed. If that is the case, the exact autoresponse 
scaling function becomes [10,11] 



and depends on the universality classes only through the 
values of the exponents a and Xr/z, while tq merely 
is a normalization constant. This prediction of local 
scale invariance (LSI) has been confirmed in a variety of 
physically very different models, see [10] and references 
therein. Still, better evidence in favour of LSI than mere 
phenomenology would be desirable. 

The origin of LSI can be understood in the special case 
z = 2 (which is realized in all cases of phase-ordering ki- 
netics we are going to consider [12]). A Ward identity can 
be written down such that if the system is known to be 
(i) scale- invariant with z = 2 and (ii) Galilei- invariant, 
then its invariance under the full group of local scale 
transformations follows [13]. Therefore, scale invariance 
and Galilei invariance appear as the building blocks for 
LSI. So far, no direct test of Galilei invariance in phase- 
ordering systems was ever performed. 

Direct tests of Galilei invariance require the study of 
the time- and space-dependence of response functions. 
For phase-ordering, it is natural to expect a scaling be- 
haviour of the linear response function 



R(t,s;r) 
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with the scaling function Fr(x, u) and Fr(x, 0) — Jr{x). 
For z = 2, Galilei invariance (combined with the Ward 
identity of LSI) predicts [14,10,13] 

/ M r 2 \ 
R(t, s; r) = R(t, s) exp — (3) 



f R (x)=r x 1+a - x */ z (x-iy 



(2) 
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where R(t, s) is the autoresponse function given by 
eqs. (1,2) and M. is a direction-dependent non-universal 
constant. Eq. (3) gives the full spatio-temporal scaling of 
the linear response. We shall sucessfully test it in the 2D 
and 31? kinetic Glauber-Ising models and provide thereby 
the first direct evidence of Galilei invariance in a phase- 
ordering system. Since the 2D/3D Glauber-Ising model 
cannot be reduced to a free-field theory, it is highly non- 
trivial that its exact response function takes the simple 
free- field form eq. (3). 

As precise simulational data for the autoresponse 
R(t, s) are difficult to obtain, it is convenient to study 
instead the thermoremanent magnetization [15] 
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TM TRM (t, s)//i (0 ) = p(t, s) = T du R(t, u) (4) 

Jo 

We consider the Ising model on a hypercubic lattice, 
with periodic boundary conditions and the Hamiltonian 
H = — <7i<7j. We use heat-bath dynamics defined 

through the stochastic rule 

<7i(t + 1) = ilwith probability \ [1 ± tanh(/i ; (t)/T)] (5) 

with the local field = X) a (i) cr a(*) an( ^ where a(i) 

runs over the nearest neighbours of the site i. -Mtrm(£, s) 
is measured by applying a quenched spatially random 
magnetic field ±/i(o) for times between the quench at 
t = and the waiting time s [15]. The presence of this 
external field then changes the local field in eq. (5) to 

M*) =Ea(i) CT aW± V)- 

It has been understood recently that the scaling be- 
haviour of Mtrm for s 3> 1 is not a simple power law 
but rather shows a cross-over behaviour [16,8]. It is of 
the form [9] 

p(t, s) = r s- a f M (t/s) + r lS - x «/ z g M {t/s) (6) 

provided the system was initially prepared at infinite 
temperature. Here the correction term can be estimated 
as 5m(^) — x~ Xr / z , while the scaling function /m(x) can 
be found from local scale invariance using eqs. (1,2,4) and 
is given explicitly by eq. (5.47) in [10]. Finally, r$ t i are 
normalization constants. In practice, it turns out that 
Xr is quite close to its lower bound d/2 and in particu- 
lar in 2D, the correction to the leading scaling behaviour 
is sizeable. Before any meaningful study of the spatio- 
temporal response can be carried out, the correction term 
must be subtracted off. 

As explained in [9], this can be carried out by fixing 
a value of x = t/s. Then r 4 are obtained by fitting 
eq. (6) to the computed thermoremancnt magnetization 
-Mtrm- In order to illustrate the quality of this proce- 
dure, we show in figure 1 data for the scaling function 
Jm{x) obtained after subtraction of the correction term 
r\s~ XR / z gM{t/ s) from the integrated autoresponse p of 
the 2D and 3D Glauber-Ising model. Statistical error 
bars in 2D are smaller than the symbol size and in 3D 
are of the order of the scatter in the data. In table I, we 
list the values of A_r [4,17] and the constants ro,i which 
were determined at x = 7. We find a nice data collapse 
and a clear agreement with the LSI prediction. 

Having thus checked the expected scaling of the au- 
toresponse function, we can now turn towards the spatio- 
temporal response. Again, we consider the integrated re- 
sponse rather than R(t, s; r) because it is considerably 
less affected by noise. In order to fix the non-universal 
parameter M. in (3), we form 

po(t,s;r) =T f dr R(t, s - r;r) (7) 
Jo 



which is measured on the lattice by keeping a small ran- 
dom magnetic field ±/i(o) until the waiting time s at the 
site and then observing the thermoremancnt magne- 
tization M TRM {t,s;r) = h( )Po(t,s;r)/T = (h <r r (t))/T 
at a different site r (where (• • •) indicates an average over 
the thermal noise whereas the bar means an average over 
the random field distribution). We computed the spa- 
tially and temporally integrated response for the Ising 
model with Glauber dynamics in two and three dimen- 
sions. The two-dimensional systems usually contained 
300 2 spins on a square lattice, whereas in three dimen- 
sions cubic lattices with typically 60 3 sites were consid- 
ered. Some Monte Carlo simulations were also done for 
other system sizes in order to check against finite-size 
effects. Since p (t,s;r) is very noisy, one must average 
over a large number of runs with different realizations of 
the thermal noise and of the spatially random magnetic 
field. For every waiting time we avaraged over at least 
50000 different runs. The whole study took about 2 • 10 5 
CPU hours on a SGI Origin 3800 parallel computer. 

Following the same lines as in [9] , we arrive at the scal- 
ing form [18], where r varies along a fixed direction on 
the lattice 

p (M;r) ^ poo(t - s:r)+r s- a F (j^^-jj 

+r 1 s- x ^ z G (^,M 1 ^j (8) 

where the universal scaling functions read 

F (x,y)= / dv exp[-y/2(x- l + v)]h R (x,v) 
Jo 

G (x,y)~x- Xli/z e- v/2x (9) 

with hji(x,v) :— /r (t 3 ^) (1 ~ w) _1_a where M is a 
direction-dependent new parameter. In this work we 
choose T such as to keep anisotropy effects small and also 
to avoid the cross-over to critical dynamics at T = T c . 
This is illustrated in figure 2, for the 2D case at T = 1.5 
where we plot Fo as function of r 2 /s. For the displayed 
data points r varies between 1 and ^fs when going along 
the (10) direction and between \[2 and V2s when go- 
ing along the (11) direction. At the considered temper- 
ature any anisotropies are very small and the values of 
M, determined separately in the two directions, coincide, 
•M(io) = Al(n) = M.. Similar results also hold in 3D at 
T = 3. We fix M at x = 7 (see table I and inset) us- 
ing (8) and then obtain a parameter-free prediction for 
other values of x. This is shown in the inset of figure 2, 
where the LSI-prediction is compared to numerical data 
for two additional values of x and two different waiting 
times. Since by now all non-universal parameters are 
fixed, this constitutes by itself a quantitative test of local 
scale invariance. 
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A fuller test of the spatio-temporal response is ob- 
tained by considering the spatially and temporally in- 
tegrated response function 



T I du / drr d_1 i2(t,ti;r) (10) 
'a J a 



where the space integral is along a straight line of length 
A = Y / 7ts but we do allow for the possibility of an 
anisotropy as a function of the solid angle Q. Such 
anisotropics are known to occur if T < T c [19]. As before, 
we derive from (3) the scaling form [18] 



dn 

+ r s a 



dp(t,s;„) =p(1){t _ sy) 

" 2 -y 2 > (t/s, n) + ri s d l 2 - XR ' z pW (t/s, n) 



(ii) 



with the explicit scaling functions 



P [2 \x,n) = 



p^(x,p) 



d Jo 

d/2 



/ dvh^Xjv)^ 
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x — 1 + v 



(12) 



This is the general expression for the scaling of the spatio- 
tcmporally integrated response function. If we fix p and 
let x = t/s vary, the form of the scaling function of 
s a ~ d / 2 dp/dQ. merely depends on p. We stress that the 
exponent \r and the free parameters r ,ri,M are now 
all fixed such that there remains no free fitting parameter 
at all when comparing (11) with our numerical data. 

It is of interest to compare the maximal distance r max 
accessed by our simulation with the physical length scale 
of the problem, for example the domain size L(s) at the 
instant the magnetic field is turned off. From the cor- 
relation function, in 2D we estimate L(s) w 3,6,9.5 for 
s = 25, 100, 225, respectively, with /i (0) = 0.05. In 3D, 
we find L(s) w 2.1, 3.3, 4 for s = 25, 64, 100. As expected 
L(s) ~ s 1 / 2 , with z « 2. The value p — 4 (p = 2) corre- 
sponds to r max /Z/(s) i=a 3.2 (3.4) in 2D (3D), respectively 
and we see that the data probe the large-distance region, 
well beyond the domain size L(s). 

The integrated spatio-temporal response in 2D at T = 
1.5 is shown in figure 3 for two values of p and where the 
correction p^ was already subtracted off. We find a nice 
scaling behaviour over the whole range of waiting times 
we could consider. We stress that the agreement between 
local scale invariance (3) and our data for several values 
of p is a real test, since no free parameter remains. In 
particular, both the height and the position of the maxi- 
mum of the scaling function for p = 2 and 4 is completely 
fixed. We can conclude that LSI is fully vindicated. This 
is the first time that the exact functional form of the scal- 
ing function of the spatio-temporal response of a generic 
non-equilibrium spin system is found. 



Similarly, the integrated spatio-temporal response in 
3D at T — 3 is shown in figure 4. Again, a nice scaling 
behaviour in perfect agreement with LSI is found. 

The prediction (3) is also verified in a few exactly solv- 
able models with a non-conserved order parameter and 
undergoing phase-ordering kinetics with z = 2. In par- 
ticular, (3) holds in the e?-dimcnsional kinetic spherical 
model with d > 2 quenched to either T < T c or T = T c 
and independently of initial correlations [20-22,16,6]. 

Summarizing, we presented the first quantitative study 
of the scaling of the spatio-temporal response in coarsen- 
ing systems. For the first time, we find direct evidence 
for Galilei invariance in the ageing regime of a phase- 
ordering system. This is a new dynamic symmetry whose 
presence in a non-equilibrium phase transition was not 
anticipated. Our result provides a strong indication that 
LSI is indeed a true space-time dynamical symmetry of 
statistical systems undergoing phase-ordering kinetics. 
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Franzosisches Hochschulzentrum (BFHZ) and by CINES 
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FIG. 1. Scaling of the integrated response function p for 
the Glauber-Ising model, (a) 2D at temperature T — 1.5 (b) 
3D at T = 3. The symbols correspond to different waiting 
times. The full curves are obtained by integrating (2). 
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FIG. 2. Scaling of the temporally integrated response func- 
tion pa for the 2D Glauber-Ising model at T = 1.5 and 
s = 225. Data obtained along the (10) resp. (11) direction 
are shown by open resp. filled symbols. Diamonds: x — 5, 
triangles: x — 7, circles: x = 9. Inset: Determination of the 
mass M, see main text. Some typical error bars are shown in 
the inset. 




FIG. 3. Scaling of the integrated spatio-temporal response 
function p for the 2D Glauber-Ising model at T = 1.5 and at 
(a) p — 1 (inset: p = 4) and (b) p = 2 (integrating (a) along 
the (10) direction and (b) along the (11) direction). The full 
curves are from eq. (12). Error bars are smaller than the 
symbol size. 
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FIG. 4. Scaling of the integrated spatio-temporal response 
function p for the 3-D Glauber-Ising model at T — 3 and at 
(a) p — 1 and (b) p = 2 (integrating (a) along the (100) di- 
rection and (b) along the (110) direction). The full curves are 
from eq. (12). Some typical error bars are also shown. 



TABLE I. Values of the autoresponse exponent Xr and of 
the parameters ro, r\ and M in the Glauber-Ising model in 
two dimensions at T = 1.5 and in three dimensions at T — 3. 







ro 


ri 


M 


2D 


1.26 


2.65 ±0.05 


-2.76 ±0.05 


4.08 ± 0.04 


3D 


1.60 


0.31 ±0.02 


0.61 ± 0.02 


4.22 ± 0.05 



4 



